Search for: All records

Abstract Polluted white dwarfs (WDs) with small surface convection zones deposit significant concentrations of heavy elements to the underlying radiative interior, presumably driving thermohaline convection. Current models of polluted WDs frequently fail to account for this effect, although its inclusion can increase the inferred accretion rate by orders of magnitude. It has been argued that this instability cannot be treated as a continuous mixing process and thus should not be considered in these models. In this work, we study 3D simulations of a thermohaline-unstable layer propagating into an underlying stable region, approximating the polluted WD scenario. We find that although thermohaline convection works to reduce driving gradients somewhat, the front continues to propagate and the system remains unstable. Importantly, the turbulent flux of metals broadly dominates over the diffusive flux in quantitative agreement with existing mixing prescriptions implemented in some stellar evolution models (except slightly below the boundary of the propagating front, where recent prescriptions neglect overshoot-like effects). Thus, our results broadly support polluted WD models that include thermohaline mixing in their estimates of the settling rate. 

Recent observations have found a growing number of hypervelocity stars with speeds of ≈1500 − 2500 km s⁻¹that could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We used an existing AREPOhydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the one-dimensional stellar evolution code MESA. We show that the supernova shock is not energetic enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the 10^5 − 6year lifetimes inferred for observed hypervelocity runaways. Although they experience a shock from a supernova less than ≈0.02 R_⊙away, our models do not experience significant interior heating, and all contract back to radii of around 0.01 R_⊙within about 10⁴years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution. 

Abstract Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in theMESAstellar evolution code. The model identifies regions in the star where the energy associated with a local super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with initial masses of 10−100M_⊙at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to ∼10⁻²M_⊙yr⁻¹can be driven by this mechanism, and occurs in a vertical band on the H-R diagram between $3.5\lesssim \log \left(T_{\mathrm{eff}}/\mathrm{K}\right)\lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses ≳20M_⊙from evolving to become red supergiants (RSGs), with the stars instead ending their lives as blue supergiants, and offers a possible explanation for the observed lack of RSGs in that mass regime. 

ABSTRACT We present a detailed modelling study of CD-30°11223 (CD-30), a hot subdwarf (sdB)-white dwarf (WD) binary identified as a double detonation supernova progenitor, using the open-source stellar evolution software MESA. We focus on implementing binary evolution models carefully tuned to match the observed characteristics of the system including log g and Teff. For the first time, we account for the structure of the hydrogen envelope throughout the modelling, and find that the inclusion of element diffusion is important for matching the observed radius and temperature. We investigate the two sdB mass solutions (0.47 and 0.54 M⊙) previously proposed for this system, strongly favouring the 0.47 M⊙ solution. The WD cooling age is compared against the sdB age using our models, which suggest an sdB likely older than the WD, contrary to the standard assumption for compact sdB-WD binaries. Subsequently, we propose a possible alternate formation channel for CD-30. We also perform binary evolution modelling of the system to study various aspects such as mass transfer, orbital period evolution, and luminosity evolution. Our models confirm CD-30 as a double detonation supernova progenitor, expected to explode ≈55 Myr from now. The WD accretes an ≈0.17 M⊙ thick helium shell that causes a detonation, leaving a 0.30 M⊙ sdB ejected at ≈750 km s−1. The final 15 Myr of the system are characterized by helium accretion which dominates the system luminosity, possibly resembling an AM CVn-type system. 

Abstract Regular, automated testing is a foundational principle of modern software development. Numerous widely used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates. 

ABSTRACT Blue large-amplitude pulsators (BLAPs) make up a rare class of hot pulsating stars with effective temperatures of ≈30 000 K and surface gravities of 4.0–5.0 dex (cgs). The evolutionary origin and current status of BLAPs is not well understood, largely based on a lack of spectroscopic observations and no available mass constraints. However, several theoretical models have been proposed that reproduce their observed properties, including studies that identify them as pulsating helium-core pre-white dwarfs (He-core pre-WDs). We present here follow-up high-speed photometry and phase-resolved spectroscopy of one of the original 14 BLAPs, OGLE-BLAP-009, discovered during the Optical Gravitational Lensing Experiment. We aim to explore its pulsation characteristics and determine stellar properties such as mass and radius in order to test the consistency of these results with He-core pre-WD models. Using the mean atmospheric parameters found using spectroscopy, we fit a spectral energy distribution to obtain a preliminary estimate of the radius, luminosity, and mass by making use of the Gaia parallax. We then compare the consistency of these results to He-core pre-WD models generated using Modules for Experiments in Stellar Astrophysics, with predicted pulsation periods implemented using gyre. We find that our mass constraints are in agreement with a low-mass He-core pre-WD of ≈0.30 M⊙. 

Abstract Type Ia supernovae (SNe Ia), critical for studying cosmic expansion, arise from thermonuclear explosions of white dwarfs, but their precise progenitor pathways remain unclear. Growing evidence supports the “double-degenerate scenario,” where two white dwarfs interact. The absence of nondegenerate companions capable of explaining the observed SN Ia rate, along with observations of hypervelocity white dwarfs, interpreted as surviving companions of such systems, provide compelling evidence for this scenario. Upcoming millihertz gravitational-wave observatories like the Laser Interferometer Space Antenna (LISA) are expected to detect thousands of double-degenerate systems, though the most compact known candidate SN Ia progenitors produce marginally detectable signals. Here, we report observations of ATLAS J1138-5139, a binary white dwarf system with an orbital period of just 28 minutes. Our analysis reveals a 1M_☉carbon–oxygen white dwarf accreting from a high-entropy helium-core white dwarf. Given its mass, the accreting carbon–oxygen white dwarf is poised to trigger a typical-luminosity SN Ia within a few million years, to evolve into a stably transferring AM Canum Venaticorum (or AM CVn) system, or undergo a merger into a massive white dwarf. ATLAS J1138-5139 provides a rare opportunity to calibrate binary evolution models by directly comparing observed orbital parameters and mass-transfer rates closer to merger than any known SN Ia progenitor. Its compact orbit ensures detectability by LISA, demonstrating the potential of millihertz gravitational-wave observatories to reveal a population of SN Ia progenitors on a Galactic scale, paving the way for multimessenger studies offering insights into the origins of these cosmologically significant explosions. 

Abstract We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range toZ= 10⁻⁶toZ= 0.2, which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity rangeZ= 10⁻⁴toZ= 0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions betweenX= 0 andX= 0.1 includingX= 10⁻², 10⁻³, 10⁻⁴, 10⁻⁵, 10⁻⁶that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables inMESAand find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from ≈20% to 80% across individual chemical mixtures, up to ≈8% and ≈15% at the bottom and top of the solar convection zone respectively, and ≈7% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity. 

Abstract Binary systems of a hot subdwarf B (sdB) star + a white dwarf (WD) with orbital periods less than 2–3 hr can come into contact due to gravitational waves and transfer mass from the sdB star to the WD before the sdB star ceases nuclear burning and contracts to become a WD. Motivated by the growing class of observed systems in this category, we study the phases of mass transfer in these systems. We find that because the residual outer hydrogen envelope accounts for a large fraction of an sdB star’s radius, sdB stars can spend a significant amount of time (∼tens of megayears) transferring this small amount of material at low rates (∼10 −10 –10 −9 M ⊙ yr −1 ) before transitioning to a phase where the bulk of their He transfers at much faster rates ( ≳10 −8 M ⊙ yr −1 ). These systems therefore spend a surprising amount of time with Roche-filling sdB donors at orbital periods longer than the range associated with He star models without an envelope. We predict that the envelope transfer phase should be detectable by searching for ellipsoidal modulation of Roche-filling objects with P orb = 30–100 minutes and T eff = 20,000–30,000 K, and that many (≥10) such systems may be found in the Galactic plane after accounting for reddening. We also argue that many of these systems may go through a phase of He transfer that matches the signatures of AM CVn systems, and that some AM CVn systems associated with young stellar populations likely descend from this channel. 

ABSTRACT We report the results from follow-up observations of two Roche-lobe filling hot subdwarf binaries with white dwarf companions predicted to have accretion discs. ZTF J213056.71+442046.5 (ZTF J2130) with a 39-min period and ZTF J205515.98+465106.5 (ZTF J2055) with a 56-min period were both discovered as subdwarf binaries with light curves that could only be explained well by including an accretion disc in their models. We performed a detailed high-resolution spectral analysis, using Keck/ESI to search for possible accretion features for both objects. We also employed polarimetric analysis using the Nordic Optical Telescope (NOT) for ZTF J2130. We did not find any signatures of an accretion disc in either object, and placed upper limits on the flux contribution and variation in degree of polarization due to the disc. Owing to the short 39-min period and availability of photometric data over 6 yr for ZTF J2130, we conducted an extensive O − C timing analysis in an attempt to look for orbital decay due to gravitational wave radiation. No such decay was detected conclusively, and a few more years of data paired with precise and consistent timing measurements were deemed necessary to constrain $$\dot{P}$$ observationally.